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Finite-Temperature Phase Structure of Lattice QCD with the Wilson 
Quark Action for Two and Four Flavors 
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We present further analyses of the finite-temperature phase structure of lattice QCD with the Wilson quark 
action based on spontaneous breakdown of parity-flavor symmetry. Results are reported on (i) an explicit demon- 
stration of spontaneous breakdown of parity-flavor symmetry beyond the critical line, (ii) phase structure and 
order of chiral transition for the case of Nf — 4 flavors, and (iii) approach toward the continuum limit. 



1. Introduction 

Elucidating the finite-temperature phase struc- 
ture of lattice QCD with the Wilson quark action 
is complicated by the fact that the action explic- 
itly breaks chiral symmetry Last year we re- 
ported an analysis of this problem[^ based on 
the idea of spontaneous breakdown of parity and 
flavor symmetry 1^]. The main findings, obtained 
through simulations of the two-fiavor system on 
an 8^ X 4 lattice, were that (i) the critical line, de- 
fined as the line of vanishing pion screening mass 
for finite temporal lattice sizes, does not extend to 
arbitrarily weak coupling, but turns back toward 
strong coupling at a finite value of (3, forming a 
cusp on the K) plane, and (ii) the thermal 
line of finite- temperature transition K = Kt{(3) 
runs past the tip of the cusp without crossing the 
critical line. 

We have also argued, based on analytical con- 
siderations including those of the two-dimensional 
Gross-Neveu model, that the cusp will move to- 
ward weak coupling as the temporal size Ni in- 
creases, pushing the thermal line in front and 
eventually pinching it at K = 1/8 and P = oo 
as Nt oo. In our view it is only in this limit 
that the true chiral phase transition emerges. 

While this view provides a consistent picture on 
how chiral phase transition arises for the Wilson 
quark action, we have left unanswered several im- 
portant questions to be clarified. These are (i) an 
explicit demonstration that parity-fiavor symme- 
try is spontaneously broken inside the cusp, (ii) 
how the phase diagram depends on the number 



of fiavors Nf , in particular if the chiral transition 
is of first order for Nf > 3 as predicted by the 
sigma model analysis, and (iii) how the tip of the 
cusp moves for an increasing temporal lattice size. 
In this article, we present results of our study on 
these questions carried out since last year. 

2. Evidence for spontaneous breakdown of 
parity-fiavor symmetry 

Spontaneous breakdown of parity-flavor sym- 
metry inside the cusp is signaled by a non- 
vanishing vacuum expectation value of the pion 
field. For the case of two flavors Nf = 2 which 
we analyze, one may take the pion condensate to 
point in the ra direction < -ipijr^T^ip >^ 0. The 
pion spectrum will then consist of a massive tt*^ 
and massless tt^ which are the Nambu-Goldstone 
modes of the broken symmetry. 

In order to ascertain these predictions we 
carry out hybrid Monte Carlo (HMC) simula- 
tions, adding a symmetry-breaking term SSw = 
2KH "^^ip^ij^T^ipn to the action to avoid in- 
frared divergences due to massless tt* modes. 

Runs are made on an 8'^ x 4 lattice at /? = 3.5 
with 0.21 < K < 0.28 in steps of 0.01. The range 
for the hopping parameter includes the interval 
Kc ^ 0.2267(2) < K < 0.2463(7) = K'^ previ- 
ously estimated to be the parity broken phase 
For the external field we take H = 0.2, 0.1, 0.05, 
0.02 and 0.01. For each value of K and H, typ- 
ically 20 — 50 trajectories with 0.5 time units 
are made for thermalization followed by 50 tra- 
jectories for measurements. The presence of ex- 
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Figure 1. Parity-flavor order parameter as a func- 
tion of 1 jK for various values of external field H. 
at /3 = 3.5 on an 8^ x 4 lattice for iV/ = 2. 

ternal field significantly reduces fiuctuations, and 
we find these statistics to be sufficient. Hadron 
screening masses are calculated by periodically 
doubling the lattice in the spatial directions. 

In fig. [l] we plot the parity-flavor order pa- 
rameter < ip'^j^Tsip > as a function of 1/K for 
H = 0.2 — 0.01. Vertical lines mark the position 
of the critical lines at /3 = 3.5. We clearly observe 
that the order parameter inside the cusp {i.e., be- 
tween the two critical values Kc and K'^) tends to 
a non-vanishing value as the external field H is 
reduced, while it decreases toward zero outside. 

The TT^ screening mass squared is plotted in 
fig. ^a) as a function of l/_fi' and in fig. ^(b) 
as a function of H. Inside the cusp de- 
creases toward zero, while outside it converges to- 
ward values calculated with H = O(open squares). 
For K = 0.23 and 0.24, which are inside the 
cusp, lines in fig. ^b) show quadratic fits of form 
^ A + BH + CH"^. A small non- vanishing 
value of the intercept A « 0.02 is ascribed to fi- 
nite spatial size effects. 

We conclude that the behavior of both the 
parity-flavor order parameter and strongly 
supports spontaneous breakdown of parity-fiavor 
symmetry inside the cusp. 

An interesting point to note is that m'^± con- 
verges to the limit H — from below for K > K'^, 
while the approach is from above for K < Kc. In 
other words, the parity-fiavor broken phase is en- 
larged beyond the upper part of the critical line 
in the presence of the external field. 

We also mention that extracting the 7r° mass 
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Figure 2. (a) m} ± as a function oiljK for various 
values of external field 7J at /3 = 3.5 on an 8'^ x 4 
lattice for the case of iV/ = 2. Open squares are 
results for H = 0. (b) Same as (a) plotted as a 
function of H for H < 0.06. 

inside the cusp requires computation of discon- 
nected quark loop contribution in the ir'^ prop- 
agator. This may be carried out with the tech- 
nique of wall source without gauge fixing previ- 
ously applied to the 77' propagator We leave 
this interesting problem for future work. 

3. Phase diagram for the case of Nf = 4 

Our analysis of the Nf — 4 system is made for 
the dual purpose of confirming the cusp structure 
of the critical line and examining the dependence 
of the order of chiral transition on Nf . 

HMC runs with H — are carried out on an 
8^ X 4 lattice at < /? < 4.0 and 0.19 < K < 0.30. 
Thermalization of 20 — 50 trajectories with 0.5 
time units followed by 50 — 100 trajectories for 
measurements are made at each (3 and K, chain- 
ing runs in K for each value of (3. The location 
of the critical line is estimated from results of 
calculated on a spatially doubled lattice, and the 
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Figure 3. Phase diagram for the iV/ = 4 sys- 
tem for Nt = 4. Previous resuh for iV/ = 2 is 
also shown for comparison. Sohd Hues are critical 
lines and dotted lines the line of thermal transi- 
tion. Solid squares mark points where first-order 
signals are found, while open rectangles represent 
region of smooth crossover. 

position and order of the thermal transition is 
examined through behavior of physical quantities 
and their time histories. 

In fig. H we show the phase diagram for TV/ = 4 
together with that for TV/ — 2 previously re- 
ported ||. We find a cusp structure similar to 
that of N j = 2 except for a shift of <5/3 « 1.2 
toward stronger coupling, which we qualitatively 
expect from a larger magnitude of sea quark ef- 
fects for N I = 4. We also find strong first-order 
signals across the thermal line away from the tip 
of the cusp as marked by solid squares. Surpris- 
ingly, however, the transition becomes weaker to- 
ward the tip of the cusp, apparently turning into 
a smooth crossover at /3 = 2.5 — 2.0. The location 
of the crossover is indicated by open rectangles in 
fig. I 

The smoothing of the transition close to the 
cusp is illustrated in fig. ^(a) for and in (b) for 
the Polyakov line 17. Jumps in physical quantities 
are apparent at /3 — 4.0 and 3.5. However, at 
/3 = 2.5 and 2.0, data taken for increasing and 
decreasing values of K overlap with each other 
and do not show any sign of discontinuity. 

At /3 = 2.0 runs are also made with an in- 
creased spatial size of 12'^ x 4. Results for 
plotted by open circles do not show deviation 
from those on an 8'^ x 4 lattice. Thus it is un- 




Figure 4. (a) as a function oil jK for various 
values /3 on an S'^ X 4 lattice forthe case of TV/ = 4. 
Open circles for (i — 2.0 are taken on an 12'^ x 
4 lattice, (b) Same as (a) for the real part of 
Polyakov loop i?ef2. 

likely that finite-size effects has rounded a first- 
order discontinuity into a smooth crossover on an 
8^ X 4 lattice. 

A possible interpretation of the smoothing is 
that it is caused by an explicit breaking of chiral 
symmry in the the Wilson quark action, whose 
effect is larger for stronger coupling. If one in- 
creases the temporal lattice size A^t, the cusp and 
the line of thermal transition will both move to- 
ward weaker coupling. Therefore the first-order 
transition may become extended up to and be- 
yond the tip of the cusp for sufficiently large At . 

Another possibility is that the first-order tran- 
sition we find is a lattice artifact, sharing its ori- 
gin with the sharpening of the A/ = 2 transition 
observed at /3 « 5.0§. This possibihty does not 
contradict the report of a first-order transition 
for A/ = 3§1. In fact the interval 4.0 < /3 < 4.7 
where first-order signals were observed is away 
from the cusp expected at /? « 3.0. Also we think 
that the divergence of runs at /? « 3.0 from a hot 
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Figure 5. Location of the critical line for the tem- 
poral sizes iVt = 4 (dotted lines) and 8 (solid 
lines) for Nf = 2 and 4. 

and a mixed starting configurations reported in 
ref. 1^ is an indication that the runs were made 
inside the parity-broken phase rather than a sig- 
nal of metastability. 

For either of the two possibilities above, if the 
chiral phase transition is of first order for Nf > 3 
in the continuum, it will emerge only for larger 
Nt for which the cusp and the thermal line move 
into the scaling region toward weak coupling. 

4. Approach toward the continuum limit 

Discussions of the previous section naturally 
raise the question how the cusp of the critical line 
moves when the temporal lattice size is increased. 
To examine this problem, we analyze the critical 
line for an 8^ x 8 lattice for Nf = 2 and 4 systems. 
We do not repeat a description of simulation de- 
tails since they are similar to those for an 8^^ x 4 
lattice. 

In fig. H we plot the critical line estimated from 
the pion mass. We observe that the cusp is clearly 
shifted toward weak coupling. However, the mag- 
nitude of shift is quite small both for Nf = 2 and 
4 so that the tip of cusp is still in the region of 
strong coupling for Nt = 8. Let us note that 
previous results by the QCDPAX Collaboration 
for Nf — 2 also indicate that the tip of the cusp 
is located at /3 = 4.0 - 4.2 for A^^ = 6 and at 

= 4.5 - 5.0 even for iVf = 18|]. A recent 
study on a symmetric lattice employing parity- 
flavor breaking external fields also shows that the 



cusp is located below /3 = 5.0 up to iVt = 100]. 
These results imply that a substantial increase in 
the temporal lattice size is needed before the cusp 
and the thermal line move into the scaling region 
(e.g., (3 > 5.5 for Nf = 2). 

We should emphasize that this result has a sig- 
nificant impact also for spectrum calculations at 
zero temperature. Since the location of the cusp 
is controlled by the smaller of the temporal and 
spatial lattice sizes, unless the spatial size is taken 
sufficiently large, critical hopping parameter will 
be absent, leading to systematic errors in mea- 
sured hadron masses. 

Overall our study indicates that within the pla- 
quette gauge action and the Wilson quark ac- 
tion large lattice sizes are needed for an explo- 
ration of continuum properties of the chiral phase 
transition. The associated computational diffi- 
culties point toward application of improvement 
ideas to achieve further progress in finite- 
temperature studies with the Wilson-type quark 
actions. 

We thank Y. Iwasaki, K. Kanaya and T. Yoshie 
for useful discussions. This work is supported in 
part by the Grants-in-Aid of the Ministry of Ed- 
ucation (Nos. 04NP0801, 08640349, 08640350). 

REFERENCES 

1. M. Fukugita, S. Ohta and A. Ukawa, Phys. 
Rev. Lett. 57 (1986) 1974; A. Ukawa, Nucl. 
Phys. B(Proc. Suppl.)9 (1990) 463. 

2. S. Aoki, A. Ukawa and T. Umemura, Phys. 
Rev. Lett. 76 (1996) 873; Nucl. Phys. B(Proc. 
Suppl.)47 (1996) 511. 

3. S. Aoki, Phys. Rev. D30 (1984) 2653; Phys. 
Rev. Lett. 57 (1986) 3136; Nucl. Phys. B314 
(1989) 79. 

4. N. Kuramashi et al, Phys. Rev. Lett. 72 
(1994) 3448. 

5. C. Bernard et al, Phys. Rev. D46 (1992) 
4741, 49 (1994) 3 574, 50 (1994) 33 77. 

6. Y. Iwasaki et al., hep-lat/9605030| . 

7. K. M. Bitar, |hep-lat/9602027| . 

8. Y. Iwasaki et al, Nucl. Phys. B(Proc. Suppl.) 
42 (1995) 502; ihid A7 (1996) 515. 

9. MILC Collaboration (presented by M. 
Wingate ), this volume. 



